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A theory is presented of DNA hairpins enclosed in a nanochannel. A hairpin becomes constrained 
as it approaches the wall of a channel which leads to an entropic force causing the hairpin to tighten. 
The free energy of the hairpin computed in the classical limit is significantly larger than what one 
would expect. As a result, the distance between hairpins or the global persistence length is often 
tens of fim long and may even reach mm sizes for 10 nm thin channels. The hairpin shape and 
size, and the DNA elongation are computed for nanoslits, and circular and square nanoschannels. 
A comparison with experiment is given. 
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I. INTRODUCTION 

There has been a longstanding interest in the effects of confining polymer chains in tubes and slits, for instance, 
in relation to the field of size exclusion chromatographyi2iiMi£i£iL&. Recent advances in chip lithography have led to 
the fabrication of well-defined nanochannels^iill so that issues in polymer physics may now be addressed which were 
elusive in the p as i|iiii2ii2iiiiSiiSiiLi2iiS. At the same time, stiff polymers like DNA and actin can now be well-aligned 
which is of considerable interest in the biosciences22i2I. 

In order to understand the mechanism of chain elongation, it is crucial to investigate the alignment as a function of 
the size of the channel cross section and this has been carried out recently by Reisner et al^. The channel width was 
of the order of the DNA persistence length. They focused especially on the relaxation times of monodisperse double- 
stranded DNA fluctuating in square and rectangular nanochannels which they compared with estimates inferred from 
theory^. The pictures of fluorescently labeled DNA they show are, however, quite remarkable with respect to their 
extreme degree of elongation, which has not been remarked upon before. Let us consider, in fact, a naive estimate of 
the distance g between two consecutive hairpins of the nanoconfined DNA which would be something like 

g = 7rrexp(U /k B T) (1) 

The radius of a hairpin, which is supposed to be semicircular, is r and its bending energy is Uq. The temperature 
is T and k B is Boltzmann's constant. A form similar to eq arises in the theory of nematically confined stiff 
chains^2i22i24i2£i2£i2L2& and will be rederived below. Note that the B-helix is assumed to remain intact. Thus, in the 
Reisner experiment, in the case of DNA confined within a 60 x 80 nm channel for instance, one would expect the 
hairpin to extend to its maximum diameter of 100 nm (equal to the diagonal). This leads to g ~ 1 /im from eq lfl|. 
This is the mere step length of the effectively one-dimensional random walk, yet the T2 DNA having a contour length 
of 63 /jm is almost fully elongated^! 

I argue here that there is an important effect missing in eq (I), namely entropic depletion of the DNA hairpin by 
the wall of the nanochannel. As the hairpin approaches the wall, it loses orientational and translational degrees of 
freedom. The resulting entropic repulsion forces the hairpin toward the axis of the channel so the hairpin is tightened 
up. Eq m seriously underestimates the elongation for this reason. My emphasis will be on an analysis of this 
phenomenon in the elastic or classical limit; an analytical solution of the relevant Fokker Planck equation for a stiff 
chain trapped in a pore is well-known to be difficult although numerical analyses have been performed22i22i2ii22. The 
hairpin is thus assumed to be a two-dimensional curve within a plane aligned along the long axis of the nanochannel. 
The hairpin chord, which is the vector distance between the two ends, is perpendicular to this axis but the shape 
of the curve is not semicircular and needs to be determined. The chord of length 2r has to be fitted into the cross 
section of the channel. This is a geometrical problem in two dimensions. Clearly, the confinement of the chord entails 
a loss in both orientational and translational degrees of freedom, giving rise to a free energy F con f of entropic origin. 
This decreases as the hairpin tightens up so it balances the bending energy U at an optimum radius f. Accordingly, 
the free energy of the hairpin in the classical limit is 

Fci = U(f) + F conf (f) (2) 

where a bar denotes the state of minimum free energy. 

My main objective will be to compute this quantity. It is convenient to express the scale g as 



g = orexp (F tot /k B T) 



(3) 
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Ftot = F cl (f) + H (4) 

Here, a is a constant and H is a sum of a variety of fiuctuational free energies beyond the classical limit which I will 
analyze qualitatively. 

The paper is organized as follows. In the next section, I compute the optimum shape of the hairpin and the 
minimum bending energy. In section 3, the free energy of confinement is evaluated for three types of cross section: 
circular, square and rectangular. The distance between hairpins or the global persistence length is calculated in section 
4. In the last section, a discussion will be given of the resulting free energies, hairpin sizes and DNA elongations in 
comparison with experiment. Therein, the terms neglected in the classical limit will also be remarked upon. 



II. HAIRPIN SHAPE 



The hairpin has a length / which is to be determined while its end-to-end distance is constrained to be 2r. It 
is described by the tangential unit vector ~u (s) at contour position s (see fig 1). The hairpin curve is assumed 
to be symmetric about the y axis which is parallel to the central axis of the nanochannel. The angle between 
~u (s) and the x axis is 7r — uj (s) so that the unit vector may be written in terms of the coordinate unit vectors 
~u = — e x cosw + ~e y sinaj. The boundary conditions are assumed to be ui (0) = ^7r and cu (1/2) — 0. 

I first compute the shape of the curve by minimizing the bending energy of the hairpin of length I 



1/2 
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U = Pk B Tj ds[—j (5) 



subject to the constraint 



1/2 A/2 

ds u ■ e x = — / ds cos uj — —r (6) 
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In eq (|SJ), P is the persistence length, T is the temperature, ks is Boltzmann's constant and duj/ds equals the inverse 
radius of curvature. Minimizing eq JSJ with respect to uj (s) leads to 



d 2 uj _ 7 
~d^ ~ 21 



2 9,2 SillW ( ? ) 



where the Lagrange multiplier 7 may be positive, zero or negative. Yamakawa and Stockmayer— have previously 
analyzed this equation in their theory of the excluded-volume effect of stiff chains but our treatment differs from 
theirs because the constraint is nonzero here (eq (0). Integration of eq (0) yields 

i2 (t) 2 = '-— < 8 » 

where the constant / must be positive. 

In the remainder of the analysis we need not compute uj (s) explicitly. The bending energy is as yet given in terms 
of three parameters 

U f 1T (9) 



Pk B T 21 I 2 

We ultimately wish to minimize this with respect to I subject to eq ©• But it turns out that it is possible to do this 
analytically via the dummy variable B defined below. 

Actually, there is another restriction on eq (JSJ : the chain cannot cross the wall of the channel. This is difficult to 
take into account mathematically because it is a local constraint defined at every point of the contour s. Nevertheless, 
on physical grounds, one expects deviations of reverse curvature away from a typical hairpin configuration to be 
forbidden in view of the increase in elastic energy associated with such bending. Accordingly, I assume the simpler 
constraint dcu/ds < holds. 

The right hand side of eq |JSj is now rewritten in terms of the new angular variable v = (jr — cu) /2 with boundary 
values v (0) = 7r/4 and v (1/2) = n/2. 

^ = -(/ + 7)^(1-5 sin 2 ^ (10) 
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The constant B < 1 is defined by 



B = -g- (id 



Note that B may be positive or negative though the combination / + 7 cannot be negative. Next, the objective is to 
the express the bending energy given by cq @ solely in terms of B. First, I integrate eq (|10|l over the entire contour 
from to 1/2. This leads to 

(/+ 7 )*=4£ (12) 

in terms of the elliptic integral 

T / 2 dv 



t/4 (1-5 sin 



2 „,\ 2 



(13) 



The constraint eq 10 is rewritten as 



V2 (ds\ 

cLu — cos lu (14) 
V dio J 



The derivative is now substituted from eq l|10[l and cos a; is expressed in terms of v. The scaled contour length a = l/r 
thus becomes a function explicitly dependent on B only. 



a" 1 ^ = £, (l - \b) - f3 



(15) 



where (3 is another elliptic integral 



r j2 1 

0= I dv(l-Bsm 2 v) 2 (16) 

Jtt/4 

The constants / and 7 can also be written in terms of B with the help of eqs (| 1 1 ft and 112|l . 

/ = 8£ 2 (2-B) (17) 



7 = 8B£ 2 (18) 

We now need to derive an expression for the bending energy which may be differentiated with respect to B as 
expediently as is possible. Eqs (fTTjl and (|T5|l are first substituted into eq © and then B is eliminated with the help 
of eq l(T5|) . 

Ur 8<f ( ^ 1 D D _^ 8/3C 



S ^7nT^ = — M = (19) 

PkbI a \ 2 J a 

This is still too complicated because it contains a. Hence, I eliminate a by inserting eq <|1(J|I into eq l|14f> and using 
eq Q12[l. Eq (|19|l then reduces to the form 

E(B)=4r)f3 (20) 

in terms of the integral 



^/^dv 9jv) r (21) 

(1 -Bsin 2 v) 2 



and the function 



g (v) = 2 sin 2 v — 1 



(22) 
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The advantage of writing the energy as the product rj(3 is that the numerator and denominator are identical in the 
respective integrands so that one may expect cancellation of terms after differentiation (apart from the factor g which, 
not inconveniently, increases monotonically with v). In fact, the derivative of eq l|20(l with respect to B introduces 
sin 2 v into the integrands which is then conveniently rewritten in terms of (l — Bsin 2 u). Terms of the kind r\(3 then 
cancel and we obtain 



1 dE 
2dB 



= B~ 



ir/2 

f3 I dv- 

?r/4 I 



1 — B sin 2 v 



3/2 



-ft 



(23) 



This form is reexpressed as a two-dimensional integral with the help of eqs (|13|1 . H16|l and (|21Jl . After rearrangement, 
this yields 



1 dE 
2dB 



tt/2 /-it/2 



it/4 J-w/i 



dudv- 



g (u) (sin 2 u — sin 2 
1-Bsin 2 u) 3/2 (1 - 5 sin 2 v) 



(24) 



Next, one notes that the function g (u) / (l — _Bsin 2 u) increases monotonically with u for all B, whatever the sign. 
The remaining factor in the integrand is antisymmetric in u and v. Hence, the derivative must be positive for all B. 
In conclusion, the bending energy is a minimum as B — > — oo. Asymptotically, we have from eqs (|13f) and l|16|) 



(25) 



with 



Ci = - In tan - = 0.881374 
8 



/3~C 2 (-S)= (26) 

with 



C 2 = \V2 



Eq (|15|l yields for the optimum length of the hairpin 



Cx 2 

The minimum elastic energy of the hairpin according to eq l|19|) is 



>/'■ = << = [%• -I) = 3.3082 (27) 



Pk B T 



E m = 4 (1 - CiC 2 ) = 1-5071 (28) 



I have integrated eqs 1|15H and (|19fl numerically in order to assess the nature of the elastic "well" (see table I). 
The optimum curve is very close to that of a semicircle (B — 0; a — it) .The bending energy only starts to increase 
markedly as B approaches unity. The hairpin is almost a straight line in that case, with highly rounded ends. 



III. ENTROPY OF DEPLETION 



In the Introduction it was argued that it is plausible to introduce a classical limit for the state of the hairpin which 
is aligned along the centerline of the nanochanncl. Then, wc need to address two two-dimensional problems, the shape 
of the hairpin treated above and how the chord of the hairpin can be fitted into the cross section of the channel. In 
the latter case, we may speak of an entropy of depletion since oricntational and translational degrees of the chain are 
cut off as the hairpin is formed and restricted by the channel walls. In this section I give approximate expressions for 
the depletion entropy for three characteristic types of nanochannels: circular, square and retangular. 
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A. Circular cross section 

The chord of the hairpin is a segment of length 2r enclosed within a circle of radius a (see fig 2). The translational 
degrees of freedom of the segment are assessed by how we may position point P in the circle which is at a distance t 
from the origin O. The orientational degrees of freedom are expressed by the angle 9. In this subsection, all lengths 
will be scaled by a for convenience. Note that we require t < 2r — 1 otherwise the segment does not fit into the circle. 
In addition, if t > 1 — 2r, the segment cannot rotate freely and the maximum angle 9 m is given by 

1 = 4r 2 + t 2 - Art cos 6 m (29) 

It is in this regime that an approximation to the entropy is derived. The fraction of realizable states is written as 

I c = (2n 2 y 1 4n [ dt t [ d9 (30) 

J2r-1 JO 

This is simply a sum over all degrees of freedom noting that 9 m — 9 m (t) . It is normalized by a factor 2ir 2 applicable 
to the state if the segment were to orient and translate freely within the circle. 

The bending forces of the previous section tend to open up the hairpin toward the channel wall. Therefore, let us 
assume r and t are close to unity and introduce the small quantities 6 and e : 5 = 1 — r and e — 1 — t (0 < e < 26). In 
that case, the angle 9 m is also small and we may write cos# m — 1 — 5$™- From eq (29) we have to the leading order 

el = 25-e (31) 
Performing the integrations in eq (30), we get the entropy of depletion 

S c = k B lnI c = h B ln S + k B ln(^^J (32) 

Although this expression is formally valid in the limit S <C 1, it turns out to be a very good approximation for all r, 
even when the segment is virtually free to rotate. As r = 1 — 6 — ► 0, I c ~ 1.200 is close to unity which is the exact 
value. Furthermore, I c divided by the integral computed numerically is a monotone increasing function of 5. The 
leading approximation to the integrand td m (t) in eq l|30|) , namely a9 m with 9 m given by eq l|31|) , is also very near the 
numerical value of t0 m (t) for all r > (1 — t) /2. 

B. Square cross section 

Trying to fit a line segment into a square is different from the case discussed above. If the side of the square has a 
length A and the length of the segment is restricted to the regime A < 2r < Ay/2, then the segment cannot rotate 
through 360 degrees (see fig 3). One of its ends (denoted by point P) is constrained to lie inside the region that is 
almost like a triangle in fig 3. Its boundary in terms of the Cartesian coordinates Xb and yb is given by 

x 2 +y 2 = 4r 2 (33) 

Again, I attempt to approximate the entropy of depletion by a suitable limiting form in this case, when the segment is 
almost as long as the diagonal of the square. For convenience, all lengths are scaled by A in the rest of this subsection. 
Thus the small quantities 5 S , e x and e y are introduced: 8 S = 1 — ry/2, e x = 1 — x and e y = 1 — y in terms of the 
coordinates x and y of point P. To the leading order, we then have from eq (33) the requirement 

e x + e v < 26 s (34) 

The orientational degree of freedom is expressed by the angle 9 which is bounded by the angles 9\ and #2 (see fig 3). 
They are given by 

cos 6*i = — (35) 
2r 

sin0 2 = ^ (36) 
It 
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Since the segment approaches a diagonal in size, it is convenient to introduce the small angles c^iand <p 2 defined by 
ipi = j7r — Q\ and ip 2 = 9 2 — jTT respectively. The left and right hand sides of eqs (35) and (36) are now expanded in 
all of the small variables which yields ipi = 8 a — e x and ip 2 — 5 S — e y to the leading order. Hence, the range of the 
orientational degree of freedom is 

9 2 - 6 1 = 2S S -e x - e y (37) 
Accordingly, the fraction of realizable states is, again, a sum over all degrees of freedom 



2 f /",,,„ „ . 2 



2S S r "ZS B -e x s c3 



I s = - / dxdy {6 2 - ft) ~ - / de x / ds y (2S S - e x - e v ) = — s - (38) 

./triangle J ^ J J 37T 



and the entropy of depletion is expressed as 

S s /k B = In I s = 3 In S s + In f^-\ (39) 

The number of corners is four which has to be included in the normalization (because of fluctuations, the chord of 
the hairpin can, of course, sample the entire square cross section of the nanochanncl) . It would appear that eq i|38|) 
is a very good approximation to the exact I s because eq (|38|l equals 0.8488 at S s = 1 which is quite close to unity. 
However, the numerically computed I s is more than twice that predicted by eq l|38|) at r = 0.5 so the approximation 
must be viewed as fairly accurate at best for all r < \\[2A. 

C. Slitlike cross section 

A simple approximation for the entropy of depletion can be given in the case of a nanoslit of large aspect ratio i.e. 
when the cross section is a rectangle of width A very much larger than height D (see fig 4a). We suppose D is so 
small that we may write 9 m ~ x/r for the maximum angle 9 m when the hairpin chord just touches one of the walls 
(see fig 4b). Ultimately, this implies that the hairpin must be sufficiently stiff (irP > D, say). If the chord or line 
segment is at a distance x from the wall, the contribution to the orientational entropy is proportional to ln(2x/7rr) 
owing to the angular restriction (9 m would be 7r/2 if the segment were to orient freely). Hence, the total orientational 
entropy is given by 

2 f D/2 ( D \ 

S or /k B = -r dxln (2x/irr) = In — - 1 (40) 

D Jo \KrJ 

The segment may also translate along the longitudinal axis of length A implying a fraction of realizable states equal 
to about (A — 2r) /A. The total entropy of depletion thus becomes 



Salit/kB = hi 



A-2r\ D 

7TT 



.4 



1 (41) 



This is only valid at small enough D. For this reason a smooth crossover from eq (|41|l to eq (|39|l does not exist as 
A—>D and r\/2 -> A. 

IV. GLOBAL PERSISTENCE LENGTH 

The distance between hairpins for a chain in a uniaxially ordered matrix scales as the global persistence length which 
is given in terms of a fluctuation theorem involving the chain susceptibility 22-26 27 . This is a rigorous way of deriving 
g but is only useful if one has a precise analytical theory for the segmental distribution inside the nanochanncl, which 
is lacking at present. I have yet to establish eqs (JIJ and 10) and in the classical limit the following physical argument 
is well known2&2&. The hairpins of length af may be viewed as defects in a one-dimensional chain model of contour 
length L. There are L/g defects of energy F c i = F c i (r) and concentration af/g. Therefore, the free energy of the 
chain becomes 

TP. t t / r \ p , 







f ar\ 


L 














k B T 
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where the first two terms derive from the ideal entropy of the defect gas. Minimization of the free energy with respect 
to g leads to 

g cl = afexp(F cl /k B T) (43) 

The free energy of confinement or depletion F con f = —TS has been discussed in the previous section. If we set S = 0, 
U = Uq and we assume the bend is semicircular, we regain eq l[T)l. Any approximations incurred in the previous 
analyses can be accounted for by adding an energy H to F c i . The global persistence length will now be evaluated for 
the various nanochannels. 



A. Circular nanochannel 



The bending energy of a hairpin is given by eqs l|19|l and (|28|l . and eq expresses the free energy of depletion 
for the hairpin enclosed in the circular nanochannel. Upon minimization, the total free energy of the hairpin in the 
classical limit 



F C d (r) _ E m P 
k B T r 



In 



3.22 



(44) 



leads to an optimum radius 



r= 3 



[ElP 2 



QE m Pa] 2 - E m P 



(45) 



Using eqs (143(1 - (I45II . I have tabulated r and g for DNA with a typical persistence length of 50 nm (see table II). These 
results will be discussed in the next section. 



B. Square nanochannel 



The entropy of depletion is now given by eq l|39|) so that the free energy of a hairpin in a square nanochannel is 



F*, (r) _ E m P 
k B T r 



-3 In 



A-ry/2 



Minimization of F£r with respect to r gives 



(E 2 m P 2 + 6\/2£ m Ap) 2 - E m P 



(46) 



(47) 



The optimum radius and distance between hairpins are displayed in table III for DNA which has a persistence length 
of 57.5 nm. These predictions will be compared with the experiments of Reisner et al. 



C. Nanoslit 



Eq H41|) leads to the free energy of a hairpin in a nanoslit 



Ff* (r) E m P 



In 



A-2r\ D 
A ) nr 



Minimization yields an optimum radius 



E m PA 
A + 2E m P 



(48) 



(49) 



Note that this expression does not depend on D. In this case, it is possible to present a convenient formula for the 
distance between hairpins 



iraE m Pr f2A + 2E m P 
9 = ^ exp 



D 



A 



(50) 



In a recent experimenljSi, the dimensions of the nanoslit were kept fixed but the persistence length was varied by 
changing the ionic strength. In table IV, g is presented as a function of P for a 100 x 1000 nm nanoslit. 
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V. DISCUSSION 



The first striking conclusion is that the distance between hairpins is considerably larger than the usual persistence 
length of DNA. This is true even when P is smaller than the typical width of the channel, irrespective of the type 
(see tables II, III and IV). In the larger channels, the effect of entropic depletion is clearly seen for the hairpins are 
tightened up significantly (compare r with a or i). At the other extreme, the global persistence length may reach 
the scale of a mm which means it may prove feasible to align bacterial chromosomes in nanochanncls completely. 
Entropic depletion is an important factor in the elongation. 

The large values of g are intimately related to the small size of the orientational fluctuations of the aligned (non- 
hairpinlike) sections of the chain with respect to the channel axis. This relation has been discussed at length for a 
wormlike chain in a nematic 27 . Here, the mathematics of a chain in a channel involves both translational and orien- 
tational degrees of freedom so is significantly more complicated. Nevertheless, the scaling estimate 7 - for the angular 
fluctuations, when amended by a small numerical coefficient (see eq ()■">. '.'>tl below), qualitatively corroborates the large 
values of g proposed here. A further implication is that the crossover to the scaling regime applicable in the limit 
of flexible chains^ is not trivial. The transition is delayed because g is still very large at P ~ channel width. The 
crossover must occur at g ~ P but, as I argue below, the classical limit is no longer valid then. 

It is interesting to note that the depletion effect is much stronger when the nanochannel is square than when it is 
circular (compare the coefficients 1.5 and 3 in eqs (|44() and (146(1 ). As a hairpin retracts from the edge of a square 
channel, its freedom to move increases rapidly. This is much less so in the case the channel is circular for the curvature 
is directed smoothly inward then. 

The entropic depletion by a nanoslit given by eq l|41|) is quite peculiar. It leads to an expression for the optimum 
radius (eq H49[) ~) which is independent of D. Moreover, for highly elongated slits, the radius never increases beyond 
the persistence length. It would be interesting to see how robust this classical computation is to fluctuations. A priori 
one might expect a semiflexible chain confined in a very thin slit of infinite dimensions to behave essentially as a 
two-dimensional chain of persistence length 2P, but this is not what is predicted here. The precise range of validity 
of eq l(4*T)l needs to be established. 

Reisner et al have presented fluorescent images of DNA in square (or almost square) nanochannelsi 3 ". The elongated 
length i? cxp of the DNA represents the typical average span of the chain and is shown in table V as a function of 
the nanochannel size in the case of T2 DNA. The span of a polymer chain is the smallest size of the box it can be 
enclosed in and in one dimension its average has been computed in the flexible limit (L 3> g\ see refs l34l and l35|) 

S = 4 9 4^ (51) 



The average span is not known for the wormlike chain model but as g tends to L or becomes larger, the root-mean- 
square extension R e should be close to S. 



R. 



2Lg-2g 2 (l-e- L lv) (52) 



Besides the hairpins, orientational fluctuations also shorten the chain. If rj is the angle between its tangent and the 
axis of the channel, we have the following average 



(r? 2 ) = 0.340 I p I (53) 

This is derived from Burkhardt's expression 3 ^ for the distribution of a harmonically confined worm, suitably (though 
approximately) rescaled to account for the hard wall repulsion in a square nanochannel (see ref l2ll) . In table V, 
the effective span (S (cos rf) or R e (cos 7/)) is compared with the experimental R ex] ji^ where the cosine is computed 
to second order via eq (|53|) . The agreement is quite good if one bears in mind that no adjustable parameters have 
been used. DNA of comparable size is also highly elongated in nanoslits of dimensions presented in table IV. This 
essentially agrees with the large values of g computed from eq l|50|) . We are currently investigating whether this 
expression is quantitatively accurate^!. 

I finally discuss a variety of fluctuational terms that have been neglected in the classical limit. 

1) Eq l|44|l may be expanded to second order in order to evaluate fluctuations in the chord or size of a hairpin. The 
root-mean-square fluctuation is then (Sr 2 ) 1 ^ 2 ~ r 2 /P. 

2) The shape of a hairpin fluctuates at constant r . The length of the hairpin is largest when the bending energy 
is a minimum. Away from the minimum, the energy and length are shown in table I from which one estimates the 
fluctuation in I to be \8l\ ~ r 2 /P. 
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3) The hairpin may tilt away from the channel axis but one expects the tilt to be given by eq (|53|l or a similar 
expression when the cross section is circular. Furthermore, the hairpin is only slightly more eccentric than semicircular, 
at best, so the influence of tilting would appear to be negligible. 

4) The entropy of depletion has been computed only approximately when the cross section of the nanochannel 
is square. In the case of DNA, this may suffice because any correction terms to eq (|39|) will be sensitive to the 
fluctuations in r discussed above. However a complete analysis may turn out to be useful for very stiff polymers like 
actin. 

5) The hairpin curve is perturbed by thermal motion. The deformations are undulatory in nature and may be 
expressed by A ~ r 3/2 /P^ 27 i 37 i 38 . 

On the whole, we conclude that the classical limit breaks down as r approaches P. For the widest channels in 
tables II and III, the predictions for g must be considered to be tentative only, for this reason. Complete analyses of 
the global persistence length including fluctuations have been given for chains in the nematic stat o 26 i 2 g . It would be 
important to have similar rigorous treatments in the case of nanochannels. 
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Tables 



B 


a 


E 


0.99999 


2.0820 


16.7402 


0.99 


2.5101 


2.7350 


0.8 


2.9035 


1.8418 


0.5 


3.0491 


1.6250 


0.1 


3.1290 


1.5772 





3.1416 


1.5708 


-0,5 


3.1846 


1.5508 


-5 


3.2705 


1.5185 


-100 


3.3058 


1.5078 


-10000 


3.3082 


1.5071 



Table I. Scaled length and scaled energy E of the hairpin as a function of B. At B = the hairpin is a semicircle 
so a — ir and E = tt/2. 



a (nm) 


10 


15 


20 


25 


30 


50 


70 


100 


150 


r (nm) 


8.6 


12.1 


15.3 


18.3 


21.1 


31 


39 


50 


65 


9 W 


2900 


199 


51 


22 


13 


4.1 


2.5 


1.8 


1.3 



Table II. Hairpin radius r and global persistence length g for DNA confined in a circular nanochannel of radius a. 
The DNA persistence length is 50 nm. 



A (nm) 


35 


70 


80 


135 


190 


370 


440 


r (nm) 


15.9 


26.0 


28.5 


40 


50 


74 


82 


g (/im) 


317 


27 


19 


7.0 


4.4 


2.5 


2.3 



Table III. Hairpin radius r and global persistence length g for DNA confined in a square nanochannel with side 
equal to A. The DNA persistence length is 57.5 nm. These results pertain to the experiments of Reisner et ali&. 



P (nm) 


(/im) 


100 


18 


150 


42 


200 


79 


250 


32 


300 


204 



Table IV. Global persistence g of DNA enclosed in a nanoslit as a function of the persistence length P. The width 
A of the nanoslit is 1000 nm and its height D is 100 nm (values pertaining to the experiments of Jo et al2i). 
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A (ran) 


35 


70 


80 


135 


190 


370 


440 


S (^tm) 








47 


38 


28 


27 


S(cos9) (/itm) 








33 


23 


12 


9 


i? e (/im) 


61 


46 


41 










R e (cos8) (/xm) 


54 


37 


33 












47 


42 


35 


23 


13 


10 


8 



Table V. Experimental elongations R exp compared with theoretical predictions R e (cos9) and S(cos6) in their 
respective regimes of validity. R eX p was measured directly from Fig 4a in refill- The T2 DNA has a contour length 
of 63 fim. 

Figures 




Fig. 1 



FIG. 1: Half of a hairpin curve in a two-dimensional Cartesian coordinate system. 
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FIG. 3: Line segment enclosed in a square cross section whose two sides represent the axes of a Cartesian Coordinate system. 
Point P is restricted to the triangular region. 
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Fig. 4 



FIG. 4: a) Line segment enclosed in the cross section of a narrow slit, b) One end of the segment touches the lower wall 
whereas its center of mass is fixed at a distance x from the wall. 



